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ABSTRACT 

Recent research in multi-valued logic for quantum computing has shown practical advantages for scaling up a 
quantum computer. [] , ] 2] Multivalued quantum systems have also been used in the framework of quantum cryp- 
tography, [4] and the concept of a qudit cluster state has been proposed by generalizing the qubit cluster state. [5] 
An evolutionary algorithm based synthesizer for ternary quantum circuits has recently been presented, [2] as well 
as a synthesis method based on matrix factorization [-3]. In this paper, a recursive synthesis method for ternary 
quantum circuits based on the Cosine-Sine unitary matrix decomposition is presented. 
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1. INTRODUCTION 

A collection of controlled three state quantum systems (qutrits) perturbed by a classical force can result in the 
state of one system controlling the evolution of a second; this is called a quantum circuit. Qutrits replace classical 
ternary bits as information units in ternary quantum computing. They are represented as a unit vector in state 
space, which is a complex three dimensional vector space, Hs- In the computational basis, the basis vectors (or 
basis states) of H3 are written in Dirac notation as |0), |1), and |2), where |0} = (1, 0, 0)-^ , |1) = (0, 1, 0)-^, and 
|1) — (0,0, 1)-^. An arbitrary vector mTis can be expressed as a linear combination \^) = uq \0)+ai |l)+a2 |2), 
ao,ai,a2 G C and |ao|^ + \ai\^ + |a2|^ = 1- The real number \ai\^ is the probability that the state vector j^f) 
will be in i^^ basis state upon measurement. Note that the basis vectors in the computational basis are ordered 
by natural numbers. The state space of ternary quantum system of n > 1 qutrits is a composite complex vector 
space formed from the algebraic tensor product 7if" of component state spaces H.^. The computational basis 
for 7^1"" consists of all possible tensor products of the computational basis vectors of the component spaces; each 
vector in this basis consists of column vectors with the entry 1 in the i-th row and zeros in all others, where 
i ranges from 1 to n. [17] An arbitrary vector I^I*) in Hf"' can be expressed as linear combination of the basis 
vectors with scalars £ C such that '^^^q \o,i\^ = 1. The real number ja^l^ is the probability that the state 
vector 1 4') will be in i^^ basis state upon measurement. 

The evolution of an n qutrit quantum system occurs via the action of a linear operator that changes the 
state vector via multiplication by a 3" x 3" unitary evolution matrix. From a computational point of view, 
the evolution matrices are quantum logic gates transforming the state vectors in Ti.f"'. From a quantum logic 
synthesis point of view, these gates need to be implemented by a universal set of quantum gates. It is a well- 
known estabhshed fact that sets of one- and two- qutrit (and in general, qudit) gates are universal [10,1]. Hence, 
logic synthesis requires that 3" x 3" evolution matrices be efficiently decomposed to the level of one and two 
qutrit gates. There are several methods from matrix theory, such as QR factorization, that have been utilized 
for this purpose in binary quantum logic synthesis. Another method is the Cosine-Sine Decomposition (CSD) 
of an arbitrary unitary matrix described in section 2. This method has been recently used by Mottonen [8] et. 
al, and Shende [7] et. al for binary logic synthesis. Recently. Bullock et.al have given a synthesis method for 
multi- valued quantum logic gates using a variation of the QR matrix factorization [•)]. This paper presents a 
CSD based method for ternary quantum logic synthesis. 
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2. THE COSINE-SINE DECOMPOSITION (CSD) 

The Cosine-Sine decomposition has been used recently [8, 7] in the synthesis of binary quantum gates, which 
are 2" x 2" unitary matrices for n qubit gates. When used in conjunction with focal optimization techniques, 
the CSD provides a recursive synthesis method with a lower number of elementary gates compared to other 
methods [7]. 

Cosine-Sine Decomposition: [9,7] Let the unitary matrix Wg c™^™ be partitioned in 2 x 2 block form 



as 
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with 2r < m. Then there exist r x r unitary matrices Ui,Vi, r x r real diagonal matrices C and S, and 
(m — r) X (m — r) unitary matrices U2, V2 such that 
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The matrices C and S are the so-called cosine-sine matrices and are of the form C=diag(cos6'i, cos 6*2, ■ • ■ ,cos6'r), 
<S=diag(sin 6*1 , sin 6*2, . . . , sinfl^), such that sin^ di + cos^ Oi = 1 for \ <i <r. 

CSD for Binary Quantum Logic Synthesis: In case of binary quantum logic, all matrices are even 
dimensional as powers of two. Hence, a given mxm unitary W can be always partitioned into m/2 x m/2 square 
blocks, giving m/2 x m/2 square matrices J7i, C/2, Vi, V2, C, S upon apphcation of the CSD. The decomposition 
for this case is given in equation (3), 




The CSD can be applied to the 2x2 block diagonal factors that occur at each iteration until one reaches the 
qubit level which involves only 2x2 matrices [8]. At each iteration level the block diagonal matrices in the 
decomposition arc realized as quantum multiplexers [7]. A quantum multiplexer is a gate acting on k-\-l qubits 
of which one is designated as the control qubit. If the control qubit is the highest order qubit, the multiplexer 
matrix is block diagonal. Depending on whether the control qubit carries |0) or the gate then performs either 
the top left block or the bottom right block of the [k + 1) x (fc + 1) block diagonal matrix on the remaining k 
bits. 
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Figure 1. A 3-qubit quantum multiplexer. The / represents two wires, one for each lower qubit. Depending on the value 
of the controlling qubit, F or G is applied to the lower two qubits. 



For instance, a quantum multiplexer matrix for 3 qubits will be a 2 x 2 block diagonal matrix with each block 
matrix of size 4x4 given in equation (4). The value of the first qubit is |0) in the location of the block matrix 
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F and |1) in the location of the block matrix G. Therefore, depending on whether the control bit carries |0} or 
|1), the gate then performs either F or G on the remaining 2 qubits respectively. 



(4) 



The cosine-sine matrices in the CSD are realized as uniformly k-controlled Ry rotations [7,8]. Such gates 
operate on 1 qubits, of which the lower k are controls and the top one is the target. A different Ry is applied 
to the target for each control bit-string. The circuit for a uniformly 2-controlled Ry rotation is given in figure 2. 
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Figure 2. A uniformly 2-controlIed Ry rotation: the lower two bits are the control bits, and the top bit is the target 
bit. In general, it requires 2*^ one qubit controlled gates to implement a uniformly fc-controUed rotation. The open circles 
represent the value |0) and the closed circles the value 

For three qubits \a) , |&), and |c), the matrix formulation of a uniformly 2-controllcd Ry rotation gate is given 
in equation (5). 



cos Oi — sin 9i 
sin 9i cos 9i 




(5) 



As |6) and |c) in equation (5) take on the values |0} and |1) in four possible combinations, Oi takes on the values 



from the set {6*0, 
arbitrary angle. 



^ij^2,^3}, resulting in different Ry gates being applied to the top most qubit. Each Oi is an 



3. CSD FOR TERNARY QUANTUM LOGIC SYNTHESIS 

For ternary quantum logic synthesis, an n-qutrit gate will be a unitary matrix W of size 3" x 3". Partition W as 
in (1) with m = 3" and r = 3"-^ so that m - r = 3" - 3""^ = 3"-^(3 - 1) = 3""^ • 2. After the application of 
the CSD, W will take the form in equation (2). The matrix blocks U2, V2 will be of size 3"^^ • 2 x 3"^^ ■ 2; hence 
an application of the CSD only on these two blocks will decompose each block into the form in equation (3). 
After these two application of the CSD and some matrix factoring, W will take the form 
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Each block matrix in the decomposition given in equations (6) - (8) above is of size 3"^^ x 3"^^. We reahze 
each block diagonal matrix as a ternary quantum multiplexer acting on n qutrits of which the highest order 
qutrit is designated as the control qutrit. Depending on which of the values |0), |1), or |2) the control qutrit 
carries, the gate then performs either the top left block, the middle block, or the bottom right block respectively 
on the remaining n — 1 qutrits. The cosine-sine matrices with identity in top-left/bottom-right block corner are 
realized as uniformly (n— l)-controllcd R^/Rz rotations. These matrices can be realized as R^ or R^ rotation [6] 
matrices in R'^ applied to the top most qutrit, controlled by the lower qutrits as they range over {|0) , |1) , |2)}. 
Each configuration of the lower qutrits leads to a different Rx or Rz gate. 

EXAMPLE Consider two qutrits being acted upon by an arbitrary gate Q. The CSD synthesis of Q is 
given in figure 3. 
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Figure 3. The decomposition of an arbitrary 2-qutrit gate Q using the CSD. Each Mi is a ternary quantum multiplexer. 
The gates {CS)x and (C5')z are uniformly l-controUed rotations 



For 1 < « < 4, each Mi gate in figure 3 is a quantum multiplexer controlled by the top qutrit and can be 
decomposed to the level of elementary gates as shown in figure 4. 
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Figure 4. Quantum Ternary Multiplexer for second qutrit and its realization in terms of Muthukrishan-Stroud gates. 
The gates labled +1 and +2 are bit shifts increasing the value of the bit by 1 and 2 mod 3 respectively. Depending on 
the value of the top control qutrit a, one of Xai is applied to the second qutrit, for a € {0, 1, 2}. 



For two qutrits, the matrix for a ternary quantum multiplexer will be a 3 x 3 block diagonal matrix given in 
equation (9). The value of the first qutrit is |0) in the location of the block matrix F, |1) in the location of the 
block matrix G, and |2) in the location of the block matrix H. Therefore, depending on whether the control bit 
carries |0), or |2), the gate then performs either F, G, or H on the remaining qutrit respectively. All blocks 
in the matrix in equation (9)are of size 3x3. 
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The gates {CS)x and {CS)z in figure 3 are uniformly l-controUcd R^ and Rz rotations respectively In either 
case, the top qutrit is controlled by the lower one, as shown in figure 5. The gate {CS)z corresponds to the 
middle matrix in equation (6). The matrix formulation of this gate as a uniformly I-controUed Rz rotation is 
given in equation (10). 
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Figure 5. Uniformly l-controUed rotations for 2-qutrits, realized as multiplexers via Muthukrishan- Stroud gates, 
gates labeled +1 and +2 are bit shifts modulo 3. 
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(10) 



Depending on the three possible binary configurations of \b), 9i takes on the values from the set {^0,^1,^2}, 
resulting in different Rz gates being applied to the top qutrit. If the lower qutrit is |0), equation (10) reduces to 
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If the lower qutrit is |1) or |2), equation (10) reduces to equations (12) and (13) respectively. 
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Hence, a 2-qutrit quantum gate can be synthesized via four 1 qutrit quantum multiplexers and three 1-qutrit 
uniformly controlled rotations on the first qutrit. In general, an n-qutrit quantum gate can be synthesized via 
four n — 1 qutrit quantum multiplexers and three uniformly n — \ controlled rotations on the top qutrit. 



4. CONCLUSIONS AND FUTURE WORK 

We give a recursive procedure for ternary quantum logic synthesis by realizing n qutrit logic gates as 3" x 3" 
unitary matrices and applying the Cosine-Sine Decomposition. We conclude that this method can synthesize a 
n qutrit gate with four multiplexers acting on n — 1 qutrits and three uniformly n — l-controUed rotations. A 
two qutrit example is given. It is our future goal to do a gate count by investigating local optimizations at each 
level of recursion. We also intend to write a CAD tool for this decomposition and get a gate count for a higher 
number of qutrits, and extend the decomposition to odd radix multi-valued quantum logic synthesis. 
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